Update of explore-JULES-ES-1p0.Rmd that uses the new packages and functions
We build and test gaussian process emulators, perform a sensitivity analyis, and constrain the input space of Jules.
Andy thinks there might be mileage in analysing the atmospheric growth, which is here. /home/h01/hadaw/Research/ES_PPE/Oct20
At the moment, this vignette is hampered by the fact that emulators are failing on a few of the outputs which represent change over the historical period. The emulator is fine for predicting absolute values in the modern period.
Andy would like to see timeseries of: cVeg, cSoil and nbp, npp in GtC and GtC/yr.
Load libraries, functions and data.
Section 2.5 in Friedlingstein et al. describes how the land carbon sink is estimated.
For each ensemble member and each variable, calculate the change from the 20 years at the start of the run, to the twenty years at the end of the run.
Initial clean of data set, removing variables that don’t change, and removing NAs (models that didn’t run).
[1] TRUE
(“probability of run failure” Could be an interesting project - logit transformation for probability of run failure? Some other ML like SVM?
There are clear run failure thresholds in the parameters rootd_ft_io and lai_max, and quite strong visual indications that a_wl_io and bio_hum_cn matter.
NEXT, histograms of failure with the LHS value
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
First, use mean NPP as an example. How does NPP respond to each parameter? NAs are removed, but zero values are still included.
Some outputs (e.g fLuc, fHarvest) have an almost perfect 1:1 relationship between modern value and change, some (nbp, npp, treeFrac) quite or moderately strong, and some (csoil, cveg) very weak or non-existant.
It appears that this ensemble is less “clear cut” in having an output that clearly distinguishes between “failed” (or close to it), and “not failed”.
Having said that, having an F0 over a threshold seems to kill the carbon cycle, as before. Here, we’ve set a threshold of 0.9 (on the normalised scale) for F0, and we remove members of the ensemble with a larger F0 than that when we build emulators.
The level 1 constraint removes any input with F0 greater than 0.9 (normalised), which removes many of the zero-carbon-cycle members up front. There are 424 ensmble members remaining.
This does constraint sequentially, which may not be a good idea.
[1] 2.410932e-09
# Plot the regular km emulator. Doesn’t look great.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 2 2 2 2 1.998154 1.997409 2 2 2 2 2 2 2 2 2 1.999201 2 1.996693 2 2 2 2 2 2 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -1842.508
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1842.5 |proj g|= 1.8652
At iterate 1 f = 1841.3 |proj g|= 1.4365
At iterate 2 f = 1839.3 |proj g|= 1.8697
At iterate 3 f = 1839 |proj g|= 1.8947
At iterate 4 f = 1839 |proj g|= 1.2143
At iterate 5 f = 1838.6 |proj g|= 1.1891
At iterate 6 f = 1838.3 |proj g|= 1.9015
At iterate 7 f = 1837.6 |proj g|= 1.871
At iterate 8 f = 1837 |proj g|= 1.8711
At iterate 9 f = 1834.9 |proj g|= 1.7796
At iterate 10 f = 1833 |proj g|= 1.8503
At iterate 11 f = 1831.1 |proj g|= 1.8644
At iterate 12 f = 1830.9 |proj g|= 1.8531
At iterate 13 f = 1829.9 |proj g|= 1.8044
At iterate 14 f = 1828.8 |proj g|= 1.8322
At iterate 15 f = 1827.7 |proj g|= 1.8111
At iterate 16 f = 1825.8 |proj g|= 1.8566
At iterate 17 f = 1822.1 |proj g|= 1.8783
At iterate 18 f = 1817.9 |proj g|= 1.8192
At iterate 19 f = 1810 |proj g|= 1.8972
At iterate 20 f = 1807.7 |proj g|= 1.3523
At iterate 21 f = 1806.7 |proj g|= 1.3328
At iterate 22 f = 1804.9 |proj g|= 1.8829
At iterate 23 f = 1803.2 |proj g|= 1.6989
At iterate 24 f = 1802.4 |proj g|= 1.7693
At iterate 25 f = 1801.2 |proj g|= 1.8451
At iterate 26 f = 1799.8 |proj g|= 1.8602
At iterate 27 f = 1799.5 |proj g|= 1.6174
At iterate 28 f = 1799.2 |proj g|= 1.5947
At iterate 29 f = 1798.8 |proj g|= 1.5503
At iterate 30 f = 1796.6 |proj g|= 1.8546
At iterate 31 f = 1795.7 |proj g|= 1.586
At iterate 32 f = 1795.4 |proj g|= 1.583
At iterate 33 f = 1795.3 |proj g|= 1.6047
At iterate 34 f = 1793.7 |proj g|= 1.5063
At iterate 35 f = 1793.4 |proj g|= 1.4599
At iterate 36 f = 1793.2 |proj g|= 1.8832
At iterate 37 f = 1793 |proj g|= 1.87
At iterate 38 f = 1792.9 |proj g|= 1.4333
At iterate 39 f = 1792.4 |proj g|= 1.6437
At iterate 40 f = 1792.1 |proj g|= 1.3598
At iterate 41 f = 1791.8 |proj g|= 1.8392
At iterate 42 f = 1791.2 |proj g|= 1.8717
At iterate 43 f = 1790.7 |proj g|= 1.8739
At iterate 44 f = 1790.5 |proj g|= 1.0022
At iterate 45 f = 1790.4 |proj g|= 1.1185
At iterate 46 f = 1790.3 |proj g|= 1.8759
At iterate 47 f = 1790.3 |proj g|= 1.875
At iterate 48 f = 1790.2 |proj g|= 1.8727
At iterate 49 f = 1790.1 |proj g|= 1.802
At iterate 50 f = 1790.1 |proj g|= 0.81547
At iterate 51 f = 1790.1 |proj g|= 0.8373
At iterate 52 f = 1790 |proj g|= 1.8728
At iterate 53 f = 1789.8 |proj g|= 0.93276
At iterate 54 f = 1789.8 |proj g|= 1.8773
At iterate 55 f = 1789.7 |proj g|= 1.8754
At iterate 56 f = 1789.7 |proj g|= 1.5203
At iterate 57 f = 1789.6 |proj g|= 1.2658
At iterate 58 f = 1789.6 |proj g|= 0.92604
At iterate 59 f = 1789.2 |proj g|= 1.876
At iterate 60 f = 1789.2 |proj g|= 1.873
At iterate 61 f = 1789.1 |proj g|= 0.80926
At iterate 62 f = 1789.1 |proj g|= 0.63617
At iterate 63 f = 1789.1 |proj g|= 1.7888
At iterate 64 f = 1789 |proj g|= 1.873
At iterate 65 f = 1789 |proj g|= 1.4746
At iterate 66 f = 1789 |proj g|= 0.71729
At iterate 67 f = 1789 |proj g|= 0.71538
At iterate 68 f = 1789 |proj g|= 0.70747
At iterate 69 f = 1789 |proj g|= 0.63382
At iterate 70 f = 1788.9 |proj g|= 0.74475
At iterate 71 f = 1788.9 |proj g|= 1.7914
At iterate 72 f = 1788.9 |proj g|= 0.95343
At iterate 73 f = 1788.8 |proj g|= 0.85278
At iterate 74 f = 1788.8 |proj g|= 0.74093
At iterate 75 f = 1788.8 |proj g|= 0.90557
At iterate 76 f = 1788.8 |proj g|= 0.29949
At iterate 77 f = 1788.8 |proj g|= 0.28806
At iterate 78 f = 1788.8 |proj g|= 0.27628
At iterate 79 f = 1788.8 |proj g|= 0.55663
At iterate 80 f = 1788.8 |proj g|= 1.7132
At iterate 81 f = 1788.8 |proj g|= 1.8713
At iterate 82 f = 1788.8 |proj g|= 1.8716
At iterate 83 f = 1788.8 |proj g|= 1.8729
At iterate 84 f = 1788.8 |proj g|= 0.83666
At iterate 85 f = 1788.8 |proj g|= 0.30137
At iterate 86 f = 1788.7 |proj g|= 0.31882
At iterate 87 f = 1788.7 |proj g|= 0.2255
At iterate 88 f = 1788.7 |proj g|= 0.42148
At iterate 89 f = 1788.7 |proj g|= 0.59752
At iterate 90 f = 1788.7 |proj g|= 0.51451
At iterate 91 f = 1788.7 |proj g|= 1.4762
At iterate 92 f = 1788.7 |proj g|= 0.69297
At iterate 93 f = 1788.7 |proj g|= 0.53596
At iterate 94 f = 1788.7 |proj g|= 0.5359
At iterate 95 f = 1788.7 |proj g|= 0.55123
At iterate 96 f = 1788.7 |proj g|= 0.37871
At iterate 97 f = 1788.7 |proj g|= 0.31548
At iterate 98 f = 1788.7 |proj g|= 0.23863
At iterate 99 f = 1788.7 |proj g|= 0.13991
At iterate 100 f = 1788.7 |proj g|= 0.13728
At iterate 101 f = 1788.7 |proj g|= 0.14306
final value 1788.683595
stopped after 101 iterations
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 2 1.995574 2 1.998154 1.790352 2 2 2 1.998084 2 2 2 2 2 1.999201 2 1.996693 1.998361 2 2 1.995738 2 1.993455 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -6192.077
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6192.1 |proj g|= 1.7779
At iterate 1 f = 6097.2 |proj g|= 1.9547
At iterate 2 f = 6092.5 |proj g|= 1.9408
At iterate 3 f = 6091.7 |proj g|= 1.94
At iterate 4 f = 6089.5 |proj g|= 1.9104
At iterate 5 f = 6081.3 |proj g|= 1.9157
At iterate 6 f = 6077.6 |proj g|= 1.9116
At iterate 7 f = 6074.3 |proj g|= 1.907
At iterate 8 f = 6073.6 |proj g|= 1.9049
At iterate 9 f = 6071.5 |proj g|= 1.8255
At iterate 10 f = 6066 |proj g|= 1.7837
At iterate 11 f = 6065.2 |proj g|= 1.7673
At iterate 12 f = 6063.4 |proj g|= 1.752
At iterate 13 f = 6062.6 |proj g|= 1.9208
At iterate 14 f = 6062.1 |proj g|= 1.7448
At iterate 15 f = 6062 |proj g|= 1.7437
At iterate 16 f = 6061.8 |proj g|= 1.7402
At iterate 17 f = 6061 |proj g|= 1.7269
At iterate 18 f = 6060.1 |proj g|= 1.7134
ys=-1.254e-02 -gs= 8.259e-01, BFGS update SKIPPED
At iterate 19 f = 6059.6 |proj g|= 1.7043
ys=-2.746e-03 -gs= 5.566e-01, BFGS update SKIPPED
At iterate 20 f = 6059.3 |proj g|= 1.7005
At iterate 21 f = 6056.3 |proj g|= 1.9169
At iterate 22 f = 6055.7 |proj g|= 1.9087
At iterate 23 f = 6055.4 |proj g|= 1.8996
At iterate 24 f = 6055.1 |proj g|= 1.6244
At iterate 25 f = 6054.9 |proj g|= 1.8958
At iterate 26 f = 6054.6 |proj g|= 1.6078
At iterate 27 f = 6054.3 |proj g|= 1.5928
At iterate 28 f = 6053.9 |proj g|= 1.8949
At iterate 29 f = 6053.8 |proj g|= 1.8942
At iterate 30 f = 6053.7 |proj g|= 1.5209
At iterate 31 f = 6053.7 |proj g|= 1.516
At iterate 32 f = 6053.3 |proj g|= 1.8931
At iterate 33 f = 6052.9 |proj g|= 1.3806
At iterate 34 f = 6052.8 |proj g|= 1.3433
At iterate 35 f = 6052 |proj g|= 1.8927
At iterate 36 f = 6051.9 |proj g|= 1.3527
At iterate 37 f = 6051.9 |proj g|= 1.0531
At iterate 38 f = 6051.8 |proj g|= 0.96661
At iterate 39 f = 6051.8 |proj g|= 0.97303
At iterate 40 f = 6051.8 |proj g|= 0.62567
At iterate 41 f = 6051.7 |proj g|= 0.68667
At iterate 42 f = 6051.7 |proj g|= 1.0212
At iterate 43 f = 6051.7 |proj g|= 0.36663
At iterate 44 f = 6051.7 |proj g|= 0.20017
At iterate 45 f = 6051.7 |proj g|= 0.20136
At iterate 46 f = 6051.7 |proj g|= 0.1629
At iterate 47 f = 6051.7 |proj g|= 0.21978
At iterate 48 f = 6051.7 |proj g|= 0.23363
At iterate 49 f = 6051.7 |proj g|= 0.17448
At iterate 50 f = 6051.7 |proj g|= 0.091728
At iterate 51 f = 6051.7 |proj g|= 0.097416
At iterate 52 f = 6051.7 |proj g|= 0.12446
At iterate 53 f = 6051.7 |proj g|= 0.07705
At iterate 54 f = 6051.7 |proj g|= 0.11677
At iterate 55 f = 6051.7 |proj g|= 0.20318
At iterate 56 f = 6051.7 |proj g|= 0.094594
At iterate 57 f = 6051.7 |proj g|= 0.11668
At iterate 58 f = 6051.7 |proj g|= 0.03121
At iterate 59 f = 6051.7 |proj g|= 0.055344
At iterate 60 f = 6051.7 |proj g|= 0.11291
At iterate 61 f = 6051.7 |proj g|= 0.066099
At iterate 62 f = 6051.7 |proj g|= 0.078269
iterations 62
function evaluations 76
segments explored during Cauchy searches 76
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0782694
final function value 6051.72
F = 6051.72
final value 6051.720344
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 2 2 2 2 1.998154 1.997409 2 2 2 2 2 2 2 2 2 1.999201 2 1.996693 2 2 2 2 2 2 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -1844.124 -1844.169 -1844.496 -1844.519
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.1 |proj g|= 1.9783
At iterate 1 f = 1827.3 |proj g|= 1.5081
At iterate 2 f = 1810.4 |proj g|= 1.9257
At iterate 3 f = 1810.1 |proj g|= 1.2762
At iterate 4 f = 1808.6 |proj g|= 1.9155
At iterate 5 f = 1803.8 |proj g|= 1.9384
At iterate 6 f = 1802.5 |proj g|= 1.9352
At iterate 7 f = 1801.9 |proj g|= 1.9195
At iterate 8 f = 1801.3 |proj g|= 1.5673
At iterate 9 f = 1801 |proj g|= 1.5746
At iterate 10 f = 1800.3 |proj g|= 1.5734
At iterate 11 f = 1800 |proj g|= 1.5698
At iterate 12 f = 1799.2 |proj g|= 1.9053
At iterate 13 f = 1799 |proj g|= 1.9022
At iterate 14 f = 1798.9 |proj g|= 1.2078
At iterate 15 f = 1798.9 |proj g|= 1.2015
At iterate 16 f = 1798.8 |proj g|= 1.4287
At iterate 17 f = 1798.3 |proj g|= 1.9001
At iterate 18 f = 1797.6 |proj g|= 1.8921
At iterate 19 f = 1796.9 |proj g|= 1.8927
At iterate 20 f = 1796.5 |proj g|= 0.9502
At iterate 21 f = 1795.9 |proj g|= 1.2526
At iterate 22 f = 1795.6 |proj g|= 1.8954
At iterate 23 f = 1795.4 |proj g|= 1.5256
At iterate 24 f = 1794.9 |proj g|= 1.8901
At iterate 25 f = 1794.7 |proj g|= 1.5025
At iterate 26 f = 1793.9 |proj g|= 1.3057
At iterate 27 f = 1793.7 |proj g|= 1.3204
At iterate 28 f = 1793.5 |proj g|= 1.4504
At iterate 29 f = 1793.4 |proj g|= 1.2978
At iterate 30 f = 1793.3 |proj g|= 1.2918
At iterate 31 f = 1793.1 |proj g|= 1.8873
At iterate 32 f = 1792.8 |proj g|= 1.8811
At iterate 33 f = 1792.6 |proj g|= 1.8802
At iterate 34 f = 1792.2 |proj g|= 1.7755
At iterate 35 f = 1791.9 |proj g|= 1.7764
At iterate 36 f = 1790.9 |proj g|= 1.7796
At iterate 37 f = 1790.8 |proj g|= 1.8808
At iterate 38 f = 1790.7 |proj g|= 1.2525
At iterate 39 f = 1790.7 |proj g|= 0.78214
At iterate 40 f = 1790.4 |proj g|= 0.82047
At iterate 41 f = 1790.3 |proj g|= 1.0794
At iterate 42 f = 1790.2 |proj g|= 1.3994
At iterate 43 f = 1790.1 |proj g|= 1.8787
At iterate 44 f = 1789.9 |proj g|= 1.2268
At iterate 45 f = 1789.9 |proj g|= 0.90265
At iterate 46 f = 1789.8 |proj g|= 1.8777
At iterate 47 f = 1789.8 |proj g|= 1.8783
At iterate 48 f = 1789.7 |proj g|= 1.5976
At iterate 49 f = 1789.7 |proj g|= 0.96825
At iterate 50 f = 1789.7 |proj g|= 1.7855
At iterate 51 f = 1789.5 |proj g|= 1.7958
At iterate 52 f = 1789.4 |proj g|= 1.7966
At iterate 53 f = 1789.4 |proj g|= 1.7892
At iterate 54 f = 1789.2 |proj g|= 1.7826
At iterate 55 f = 1789.2 |proj g|= 1.194
At iterate 56 f = 1789.2 |proj g|= 1.1932
At iterate 57 f = 1789.1 |proj g|= 1.1212
At iterate 58 f = 1789 |proj g|= 1.1135
At iterate 59 f = 1789 |proj g|= 1.875
At iterate 60 f = 1789 |proj g|= 1.7381
At iterate 61 f = 1788.9 |proj g|= 1.4756
At iterate 62 f = 1788.9 |proj g|= 1.7829
At iterate 63 f = 1788.8 |proj g|= 1.7807
At iterate 64 f = 1788.8 |proj g|= 1.7842
At iterate 65 f = 1788.8 |proj g|= 0.55799
At iterate 66 f = 1788.8 |proj g|= 0.54893
At iterate 67 f = 1788.8 |proj g|= 1.162
At iterate 68 f = 1788.8 |proj g|= 1.8723
At iterate 69 f = 1788.8 |proj g|= 1.0113
At iterate 70 f = 1788.8 |proj g|= 0.42081
At iterate 71 f = 1788.8 |proj g|= 0.40831
At iterate 72 f = 1788.8 |proj g|= 0.37755
At iterate 73 f = 1788.7 |proj g|= 0.55045
At iterate 74 f = 1788.7 |proj g|= 1.8706
At iterate 75 f = 1788.7 |proj g|= 0.62034
At iterate 76 f = 1788.7 |proj g|= 0.32525
At iterate 77 f = 1788.7 |proj g|= 0.30756
At iterate 78 f = 1788.7 |proj g|= 0.39267
At iterate 79 f = 1788.7 |proj g|= 0.63073
At iterate 80 f = 1788.7 |proj g|= 0.54849
At iterate 81 f = 1788.7 |proj g|= 0.37813
At iterate 82 f = 1788.7 |proj g|= 0.28969
At iterate 83 f = 1788.7 |proj g|= 0.40478
At iterate 84 f = 1788.7 |proj g|= 0.5533
At iterate 85 f = 1788.7 |proj g|= 0.68968
At iterate 86 f = 1788.7 |proj g|= 1.2444
At iterate 87 f = 1788.7 |proj g|= 0.099835
At iterate 88 f = 1788.7 |proj g|= 0.099618
At iterate 89 f = 1788.7 |proj g|= 0.1305
At iterate 90 f = 1788.7 |proj g|= 0.13061
At iterate 91 f = 1788.7 |proj g|= 0.47715
At iterate 92 f = 1788.7 |proj g|= 0.21897
At iterate 93 f = 1788.7 |proj g|= 0.1916
At iterate 94 f = 1788.7 |proj g|= 0.22821
At iterate 95 f = 1788.7 |proj g|= 1.0767
At iterate 96 f = 1788.7 |proj g|= 0.3527
At iterate 97 f = 1788.7 |proj g|= 0.33189
At iterate 98 f = 1788.7 |proj g|= 0.21486
At iterate 99 f = 1788.7 |proj g|= 0.1876
At iterate 100 f = 1788.7 |proj g|= 0.22777
At iterate 101 f = 1788.7 |proj g|= 0.53975
final value 1788.683283
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.2 |proj g|= 1.901
At iterate 1 f = 1829.2 |proj g|= 1.7095
At iterate 2 f = 1826.3 |proj g|= 1.7222
At iterate 3 f = 1817.4 |proj g|= 1.8951
At iterate 4 f = 1810.3 |proj g|= 1.9269
At iterate 5 f = 1809.7 |proj g|= 1.3481
At iterate 6 f = 1808.7 |proj g|= 1.3548
At iterate 7 f = 1807.2 |proj g|= 1.9013
At iterate 8 f = 1805.9 |proj g|= 1.9149
At iterate 9 f = 1805.4 |proj g|= 1.8995
At iterate 10 f = 1803.8 |proj g|= 1.1077
At iterate 11 f = 1803.5 |proj g|= 1.9043
At iterate 12 f = 1803.4 |proj g|= 1.0665
At iterate 13 f = 1802.9 |proj g|= 1.9108
At iterate 14 f = 1802.7 |proj g|= 1.0578
At iterate 15 f = 1801.1 |proj g|= 1.9175
At iterate 16 f = 1800.8 |proj g|= 1.6693
At iterate 17 f = 1800.1 |proj g|= 1.6028
At iterate 18 f = 1799.4 |proj g|= 1.9058
At iterate 19 f = 1798 |proj g|= 1.8989
At iterate 20 f = 1797.2 |proj g|= 0.94213
At iterate 21 f = 1796.9 |proj g|= 1.8982
At iterate 22 f = 1796.3 |proj g|= 1.3185
At iterate 23 f = 1794.8 |proj g|= 1.524
At iterate 24 f = 1794.4 |proj g|= 1.5174
At iterate 25 f = 1793.8 |proj g|= 1.8934
At iterate 26 f = 1793.7 |proj g|= 1.8903
At iterate 27 f = 1793.6 |proj g|= 1.4863
At iterate 28 f = 1793.4 |proj g|= 1.4153
At iterate 29 f = 1793.3 |proj g|= 1.6275
At iterate 30 f = 1793.1 |proj g|= 1.7603
At iterate 31 f = 1792.9 |proj g|= 1.8829
At iterate 32 f = 1792.7 |proj g|= 1.8851
At iterate 33 f = 1792.5 |proj g|= 1.7591
At iterate 34 f = 1792.4 |proj g|= 0.74821
At iterate 35 f = 1792 |proj g|= 0.67033
At iterate 36 f = 1791.2 |proj g|= 0.88902
At iterate 37 f = 1791 |proj g|= 1.8859
At iterate 38 f = 1790.9 |proj g|= 0.8918
At iterate 39 f = 1790.9 |proj g|= 1.885
At iterate 40 f = 1790.9 |proj g|= 0.89031
At iterate 41 f = 1790.8 |proj g|= 0.88962
At iterate 42 f = 1790.8 |proj g|= 1.884
At iterate 43 f = 1790.7 |proj g|= 1.6486
At iterate 44 f = 1790.7 |proj g|= 0.87749
At iterate 45 f = 1790.5 |proj g|= 0.87374
At iterate 46 f = 1790.2 |proj g|= 0.92403
At iterate 47 f = 1790 |proj g|= 0.89554
At iterate 48 f = 1789.9 |proj g|= 0.85871
At iterate 49 f = 1789.8 |proj g|= 0.81672
At iterate 50 f = 1789.8 |proj g|= 1.7848
At iterate 51 f = 1789.8 |proj g|= 1.7893
At iterate 52 f = 1789.7 |proj g|= 1.7884
At iterate 53 f = 1789.7 |proj g|= 0.75862
At iterate 54 f = 1789.7 |proj g|= 0.7675
At iterate 55 f = 1789.6 |proj g|= 0.77272
At iterate 56 f = 1789.5 |proj g|= 1.8809
At iterate 57 f = 1789.4 |proj g|= 0.51739
At iterate 58 f = 1789.3 |proj g|= 1.8875
At iterate 59 f = 1789.2 |proj g|= 1.8826
At iterate 60 f = 1789.2 |proj g|= 1.7779
At iterate 61 f = 1789.2 |proj g|= 1.3471
At iterate 62 f = 1789.1 |proj g|= 0.54392
At iterate 63 f = 1789.1 |proj g|= 0.86088
At iterate 64 f = 1789.1 |proj g|= 1.8811
At iterate 65 f = 1789.1 |proj g|= 1.8811
At iterate 66 f = 1789 |proj g|= 1.8798
At iterate 67 f = 1788.9 |proj g|= 1.8743
At iterate 68 f = 1788.9 |proj g|= 0.80123
At iterate 69 f = 1788.9 |proj g|= 0.78927
At iterate 70 f = 1788.8 |proj g|= 1.8759
At iterate 71 f = 1788.8 |proj g|= 0.71153
At iterate 72 f = 1788.7 |proj g|= 0.52376
At iterate 73 f = 1788.7 |proj g|= 0.37546
At iterate 74 f = 1788.7 |proj g|= 0.32805
At iterate 75 f = 1788.7 |proj g|= 0.31573
At iterate 76 f = 1788.7 |proj g|= 0.29298
At iterate 77 f = 1788.7 |proj g|= 1.8724
At iterate 78 f = 1788.7 |proj g|= 0.77072
At iterate 79 f = 1788.7 |proj g|= 0.42465
At iterate 80 f = 1788.7 |proj g|= 0.30258
At iterate 81 f = 1788.7 |proj g|= 0.55782
At iterate 82 f = 1788.7 |proj g|= 0.59316
At iterate 83 f = 1788.7 |proj g|= 0.61489
At iterate 84 f = 1788.7 |proj g|= 0.50715
At iterate 85 f = 1788.7 |proj g|= 0.16702
At iterate 86 f = 1788.7 |proj g|= 0.28302
At iterate 87 f = 1788.7 |proj g|= 0.33463
At iterate 88 f = 1788.7 |proj g|= 0.67044
At iterate 89 f = 1788.7 |proj g|= 0.66924
At iterate 90 f = 1788.7 |proj g|= 0.21373
At iterate 91 f = 1788.7 |proj g|= 0.21343
At iterate 92 f = 1788.7 |proj g|= 0.12906
At iterate 93 f = 1788.7 |proj g|= 0.23661
At iterate 94 f = 1788.7 |proj g|= 0.15096
At iterate 95 f = 1788.7 |proj g|= 0.074849
At iterate 96 f = 1788.7 |proj g|= 0.089774
At iterate 97 f = 1788.7 |proj g|= 0.15525
At iterate 98 f = 1788.7 |proj g|= 0.20395
At iterate 99 f = 1788.7 |proj g|= 0.19495
At iterate 100 f = 1788.7 |proj g|= 0.12942
At iterate 101 f = 1788.7 |proj g|= 0.12935
final value 1788.682106
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.5 |proj g|= 1.8965
At iterate 1 f = 1837.2 |proj g|= 1.4588
At iterate 2 f = 1830.9 |proj g|= 1.6454
ys=-1.833e-01 -gs= 5.897e+00, BFGS update SKIPPED
At iterate 3 f = 1824.7 |proj g|= 1.8975
At iterate 4 f = 1821.1 |proj g|= 1.8605
ys=-6.464e-01 -gs= 3.354e+00, BFGS update SKIPPED
At iterate 5 f = 1815 |proj g|= 1.9047
At iterate 6 f = 1810.3 |proj g|= 1.8804
At iterate 7 f = 1809.5 |proj g|= 1.5698
At iterate 8 f = 1808.6 |proj g|= 1.8781
At iterate 9 f = 1808 |proj g|= 1.8382
At iterate 10 f = 1807.5 |proj g|= 1.657
At iterate 11 f = 1806.9 |proj g|= 1.8304
At iterate 12 f = 1806.3 |proj g|= 1.8426
At iterate 13 f = 1804.9 |proj g|= 1.8294
At iterate 14 f = 1804.6 |proj g|= 1.6177
At iterate 15 f = 1804.3 |proj g|= 1.5993
At iterate 16 f = 1804.1 |proj g|= 1.853
At iterate 17 f = 1803.7 |proj g|= 1.568
At iterate 18 f = 1803.5 |proj g|= 1.8548
At iterate 19 f = 1798.5 |proj g|= 1.8698
At iterate 20 f = 1798 |proj g|= 1.5165
At iterate 21 f = 1796.5 |proj g|= 1.8676
At iterate 22 f = 1796.4 |proj g|= 1.5568
At iterate 23 f = 1796.2 |proj g|= 1.5642
At iterate 24 f = 1793.6 |proj g|= 1.5362
At iterate 25 f = 1793.3 |proj g|= 1.8045
At iterate 26 f = 1793.1 |proj g|= 1.8804
At iterate 27 f = 1792.9 |proj g|= 1.8841
At iterate 28 f = 1792.6 |proj g|= 1.8834
At iterate 29 f = 1792 |proj g|= 1.4467
At iterate 30 f = 1791.8 |proj g|= 1.2888
At iterate 31 f = 1790.9 |proj g|= 1.8787
At iterate 32 f = 1790.7 |proj g|= 1.0173
At iterate 33 f = 1790.5 |proj g|= 0.98947
At iterate 34 f = 1790.1 |proj g|= 1.8805
At iterate 35 f = 1790 |proj g|= 1.4888
At iterate 36 f = 1790 |proj g|= 0.70684
At iterate 37 f = 1790 |proj g|= 0.91985
At iterate 38 f = 1790 |proj g|= 1.6393
At iterate 39 f = 1789.8 |proj g|= 1.882
At iterate 40 f = 1789.7 |proj g|= 1.881
At iterate 41 f = 1789.6 |proj g|= 1.8776
At iterate 42 f = 1789.6 |proj g|= 0.76916
At iterate 43 f = 1789.5 |proj g|= 0.71971
At iterate 44 f = 1789.5 |proj g|= 1.1675
At iterate 45 f = 1789.5 |proj g|= 1.6044
At iterate 46 f = 1789.4 |proj g|= 1.8782
At iterate 47 f = 1789.2 |proj g|= 1.8816
At iterate 48 f = 1789.2 |proj g|= 1.801
At iterate 49 f = 1789.1 |proj g|= 1.7905
At iterate 50 f = 1789.1 |proj g|= 0.92407
At iterate 51 f = 1789.1 |proj g|= 0.83788
At iterate 52 f = 1789 |proj g|= 0.77994
At iterate 53 f = 1789 |proj g|= 0.60603
At iterate 54 f = 1788.9 |proj g|= 0.5705
At iterate 55 f = 1788.9 |proj g|= 1.8185
At iterate 56 f = 1788.9 |proj g|= 1.0668
At iterate 57 f = 1788.9 |proj g|= 1.0849
At iterate 58 f = 1788.9 |proj g|= 0.99779
At iterate 59 f = 1788.8 |proj g|= 0.30721
At iterate 60 f = 1788.8 |proj g|= 0.24539
At iterate 61 f = 1788.8 |proj g|= 1.0934
At iterate 62 f = 1788.8 |proj g|= 1.1857
At iterate 63 f = 1788.8 |proj g|= 0.40884
At iterate 64 f = 1788.8 |proj g|= 0.35038
At iterate 65 f = 1788.8 |proj g|= 0.57915
At iterate 66 f = 1788.8 |proj g|= 0.55331
At iterate 67 f = 1788.8 |proj g|= 0.75052
At iterate 68 f = 1788.7 |proj g|= 0.99308
At iterate 69 f = 1788.7 |proj g|= 1.7903
At iterate 70 f = 1788.7 |proj g|= 0.57204
At iterate 71 f = 1788.7 |proj g|= 0.30753
At iterate 72 f = 1788.7 |proj g|= 0.40401
At iterate 73 f = 1788.7 |proj g|= 1.4657
At iterate 74 f = 1788.7 |proj g|= 0.55761
At iterate 75 f = 1788.7 |proj g|= 0.49068
At iterate 76 f = 1788.7 |proj g|= 0.42876
At iterate 77 f = 1788.7 |proj g|= 0.46333
At iterate 78 f = 1788.7 |proj g|= 1.208
At iterate 79 f = 1788.7 |proj g|= 0.6762
At iterate 80 f = 1788.7 |proj g|= 0.49907
At iterate 81 f = 1788.7 |proj g|= 0.33404
At iterate 82 f = 1788.7 |proj g|= 0.16195
At iterate 83 f = 1788.7 |proj g|= 0.40693
At iterate 84 f = 1788.7 |proj g|= 0.19521
At iterate 85 f = 1788.7 |proj g|= 0.16413
At iterate 86 f = 1788.7 |proj g|= 0.13068
At iterate 87 f = 1788.7 |proj g|= 0.061973
At iterate 88 f = 1788.7 |proj g|= 0.1542
At iterate 89 f = 1788.7 |proj g|= 0.18834
At iterate 90 f = 1788.7 |proj g|= 0.11761
At iterate 91 f = 1788.7 |proj g|= 0.11372
At iterate 92 f = 1788.7 |proj g|= 0.13048
At iterate 93 f = 1788.7 |proj g|= 0.41079
At iterate 94 f = 1788.7 |proj g|= 0.13047
At iterate 95 f = 1788.7 |proj g|= 0.060687
At iterate 96 f = 1788.7 |proj g|= 0.061571
At iterate 97 f = 1788.7 |proj g|= 0.13531
At iterate 98 f = 1788.7 |proj g|= 0.20686
At iterate 99 f = 1788.7 |proj g|= 0.13803
At iterate 100 f = 1788.7 |proj g|= 0.44702
At iterate 101 f = 1788.7 |proj g|= 0.17965
final value 1788.680770
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.5 |proj g|= 1.9627
At iterate 1 f = 1841.8 |proj g|= 1.5053
At iterate 2 f = 1835.2 |proj g|= 1.5189
At iterate 3 f = 1832.8 |proj g|= 1.9281
At iterate 4 f = 1828 |proj g|= 1.2229
At iterate 5 f = 1823.6 |proj g|= 1.0775
At iterate 6 f = 1822 |proj g|= 1.9408
At iterate 7 f = 1809.7 |proj g|= 1.8918
At iterate 8 f = 1803 |proj g|= 1.845
At iterate 9 f = 1802.6 |proj g|= 1.7398
At iterate 10 f = 1801.4 |proj g|= 1.2403
At iterate 11 f = 1799.9 |proj g|= 1.8317
At iterate 12 f = 1799 |proj g|= 1.8231
At iterate 13 f = 1797.6 |proj g|= 0.99073
At iterate 14 f = 1795.2 |proj g|= 1.0283
At iterate 15 f = 1793.4 |proj g|= 1.8782
At iterate 16 f = 1792.8 |proj g|= 1.7883
At iterate 17 f = 1792.7 |proj g|= 1.544
At iterate 18 f = 1792.6 |proj g|= 1.8611
At iterate 19 f = 1792.5 |proj g|= 1.8626
At iterate 20 f = 1791.9 |proj g|= 1.867
At iterate 21 f = 1791.7 |proj g|= 1.4722
At iterate 22 f = 1791.6 |proj g|= 1.456
At iterate 23 f = 1791.3 |proj g|= 1.7959
At iterate 24 f = 1791 |proj g|= 1.7963
At iterate 25 f = 1790.9 |proj g|= 1.8665
At iterate 26 f = 1790.7 |proj g|= 1.8672
At iterate 27 f = 1790.7 |proj g|= 0.78433
At iterate 28 f = 1790.6 |proj g|= 0.78289
At iterate 29 f = 1790.6 |proj g|= 0.7744
At iterate 30 f = 1790.4 |proj g|= 1.8684
At iterate 31 f = 1790.4 |proj g|= 1.7928
At iterate 32 f = 1790.3 |proj g|= 0.80236
At iterate 33 f = 1790.3 |proj g|= 0.82648
At iterate 34 f = 1790.1 |proj g|= 0.87867
At iterate 35 f = 1790 |proj g|= 0.8733
At iterate 36 f = 1789.8 |proj g|= 1.7905
At iterate 37 f = 1789.7 |proj g|= 1.8743
At iterate 38 f = 1789.7 |proj g|= 1.8715
At iterate 39 f = 1789.6 |proj g|= 0.80037
At iterate 40 f = 1789.6 |proj g|= 0.80415
At iterate 41 f = 1789.5 |proj g|= 0.80189
At iterate 42 f = 1789.5 |proj g|= 0.78967
At iterate 43 f = 1789.4 |proj g|= 0.75391
At iterate 44 f = 1789.3 |proj g|= 1.872
At iterate 45 f = 1789.3 |proj g|= 0.68717
At iterate 46 f = 1789.1 |proj g|= 1.4814
At iterate 47 f = 1789.1 |proj g|= 1.487
At iterate 48 f = 1789 |proj g|= 1.8742
At iterate 49 f = 1789 |proj g|= 1.4876
At iterate 50 f = 1788.9 |proj g|= 1.7243
At iterate 51 f = 1788.9 |proj g|= 0.49094
At iterate 52 f = 1788.9 |proj g|= 0.5966
At iterate 53 f = 1788.9 |proj g|= 0.59535
At iterate 54 f = 1788.9 |proj g|= 1.8733
At iterate 55 f = 1788.8 |proj g|= 0.93604
At iterate 56 f = 1788.8 |proj g|= 0.85211
At iterate 57 f = 1788.8 |proj g|= 1.4377
At iterate 58 f = 1788.8 |proj g|= 1.4772
At iterate 59 f = 1788.8 |proj g|= 0.78361
At iterate 60 f = 1788.8 |proj g|= 1.8383
At iterate 61 f = 1788.8 |proj g|= 0.3324
At iterate 62 f = 1788.8 |proj g|= 0.29873
At iterate 63 f = 1788.8 |proj g|= 0.94247
At iterate 64 f = 1788.8 |proj g|= 0.30321
At iterate 65 f = 1788.8 |proj g|= 0.32765
At iterate 66 f = 1788.7 |proj g|= 0.48866
At iterate 67 f = 1788.7 |proj g|= 1.1055
At iterate 68 f = 1788.7 |proj g|= 0.71111
At iterate 69 f = 1788.7 |proj g|= 0.50294
At iterate 70 f = 1788.7 |proj g|= 0.95867
At iterate 71 f = 1788.7 |proj g|= 0.49225
At iterate 72 f = 1788.7 |proj g|= 0.33465
At iterate 73 f = 1788.7 |proj g|= 0.3246
At iterate 74 f = 1788.7 |proj g|= 0.27462
At iterate 75 f = 1788.7 |proj g|= 0.22952
At iterate 76 f = 1788.7 |proj g|= 0.47635
At iterate 77 f = 1788.7 |proj g|= 0.36506
At iterate 78 f = 1788.7 |proj g|= 0.2705
At iterate 79 f = 1788.7 |proj g|= 0.35958
At iterate 80 f = 1788.7 |proj g|= 0.21413
At iterate 81 f = 1788.7 |proj g|= 0.083278
At iterate 82 f = 1788.7 |proj g|= 0.069125
At iterate 83 f = 1788.7 |proj g|= 0.044831
At iterate 84 f = 1788.7 |proj g|= 0.13009
At iterate 85 f = 1788.7 |proj g|= 0.083972
At iterate 86 f = 1788.7 |proj g|= 0.038698
At iterate 87 f = 1788.7 |proj g|= 0.055152
At iterate 88 f = 1788.7 |proj g|= 0.02741
At iterate 89 f = 1788.7 |proj g|= 0.026779
At iterate 90 f = 1788.7 |proj g|= 0.0423
At iterate 91 f = 1788.7 |proj g|= 0.065137
At iterate 92 f = 1788.7 |proj g|= 0.046154
At iterate 93 f = 1788.7 |proj g|= 0.041203
At iterate 94 f = 1788.7 |proj g|= 0.058948
At iterate 95 f = 1788.7 |proj g|= 0.015991
At iterate 96 f = 1788.7 |proj g|= 0.021177
iterations 96
function evaluations 123
segments explored during Cauchy searches 101
BFGS updates skipped 0
active bounds at final generalized Cauchy point 23
norm of the final projected gradient 0.0211773
final function value 1788.68
F = 1788.68
final value 1788.679245
converged
* The 4 best values (multistart) obtained are:
1788.683 1788.682 1788.681 1788.679
* The model corresponding to the best one (1788.679) is stored.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 2 1.995574 2 1.998154 1.790352 2 2 2 1.998084 2 2 2 2 2 1.999201 2 1.996693 1.998361 2 2 1.995738 2 1.993455 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -6193.86 -6193.966 -6195.06 -6196.14
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6193.9 |proj g|= 1.9414
At iterate 1 f = 6099.8 |proj g|= 1.5695
At iterate 2 f = 6091.5 |proj g|= 1.903
At iterate 3 f = 6089.3 |proj g|= 1.8889
At iterate 4 f = 6083.7 |proj g|= 1.8194
At iterate 5 f = 6082.8 |proj g|= 1.8071
At iterate 6 f = 6079.7 |proj g|= 1.9016
At iterate 7 f = 6077.9 |proj g|= 1.8974
At iterate 8 f = 6068.1 |proj g|= 1.4921
At iterate 9 f = 6066.4 |proj g|= 1.4399
At iterate 10 f = 6064.5 |proj g|= 1.9119
At iterate 11 f = 6063.5 |proj g|= 1.5745
At iterate 12 f = 6060.4 |proj g|= 1.5219
At iterate 13 f = 6060 |proj g|= 1.5145
At iterate 14 f = 6059 |proj g|= 1.9074
At iterate 15 f = 6057.9 |proj g|= 1.8963
At iterate 16 f = 6056.2 |proj g|= 1.2182
At iterate 17 f = 6054.7 |proj g|= 1.8811
At iterate 18 f = 6053.9 |proj g|= 1.1028
At iterate 19 f = 6053.2 |proj g|= 1.0823
At iterate 20 f = 6053.1 |proj g|= 1.888
At iterate 21 f = 6052.9 |proj g|= 1.0711
At iterate 22 f = 6052.7 |proj g|= 1.4231
At iterate 23 f = 6052.5 |proj g|= 1.8888
At iterate 24 f = 6052.4 |proj g|= 1.3347
At iterate 25 f = 6052.4 |proj g|= 1.0392
At iterate 26 f = 6052.4 |proj g|= 1.8869
At iterate 27 f = 6052.3 |proj g|= 1.0399
At iterate 28 f = 6052.3 |proj g|= 1.0381
At iterate 29 f = 6052.3 |proj g|= 1.023
At iterate 30 f = 6052.1 |proj g|= 0.80119
At iterate 31 f = 6051.9 |proj g|= 0.74419
At iterate 32 f = 6051.9 |proj g|= 0.69194
At iterate 33 f = 6051.8 |proj g|= 1.7041
At iterate 34 f = 6051.8 |proj g|= 0.83632
At iterate 35 f = 6051.8 |proj g|= 0.83293
At iterate 36 f = 6051.7 |proj g|= 0.64953
At iterate 37 f = 6051.7 |proj g|= 0.48825
At iterate 38 f = 6051.7 |proj g|= 0.28695
At iterate 39 f = 6051.7 |proj g|= 1.6661
At iterate 40 f = 6051.7 |proj g|= 0.46224
At iterate 41 f = 6051.7 |proj g|= 0.15696
At iterate 42 f = 6051.7 |proj g|= 0.14672
At iterate 43 f = 6051.7 |proj g|= 0.11677
At iterate 44 f = 6051.7 |proj g|= 0.048306
At iterate 45 f = 6051.7 |proj g|= 0.048253
At iterate 46 f = 6051.7 |proj g|= 0.09721
At iterate 47 f = 6051.7 |proj g|= 0.11672
At iterate 48 f = 6051.7 |proj g|= 0.103
At iterate 49 f = 6051.7 |proj g|= 0.10676
At iterate 50 f = 6051.7 |proj g|= 0.012401
iterations 50
function evaluations 56
segments explored during Cauchy searches 62
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0124006
final function value 6051.72
F = 6051.72
final value 6051.720344
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6194 |proj g|= 1.896
At iterate 1 f = 6096.2 |proj g|= 1.9139
At iterate 2 f = 6095.7 |proj g|= 1.9117
At iterate 3 f = 6092.7 |proj g|= 1.9101
At iterate 4 f = 6092.3 |proj g|= 1.9112
At iterate 5 f = 6091.6 |proj g|= 1.9106
ys=-1.660e-03 -gs= 7.046e-01, BFGS update SKIPPED
At iterate 6 f = 6088.6 |proj g|= 1.9268
At iterate 7 f = 6080 |proj g|= 1.9345
ys=-1.615e-02 -gs= 8.492e+00, BFGS update SKIPPED
At iterate 8 f = 6077 |proj g|= 1.9229
At iterate 9 f = 6069.1 |proj g|= 1.8154
At iterate 10 f = 6068.5 |proj g|= 1.6384
At iterate 11 f = 6068 |proj g|= 1.9184
At iterate 12 f = 6067.9 |proj g|= 1.9211
At iterate 13 f = 6067 |proj g|= 1.9248
At iterate 14 f = 6063.9 |proj g|= 1.9201
At iterate 15 f = 6062.4 |proj g|= 1.9124
At iterate 16 f = 6061.1 |proj g|= 1.7501
At iterate 17 f = 6059.5 |proj g|= 1.9052
At iterate 18 f = 6059.1 |proj g|= 1.7219
At iterate 19 f = 6057.7 |proj g|= 1.7267
At iterate 20 f = 6056.9 |proj g|= 1.7136
At iterate 21 f = 6056.1 |proj g|= 1.6472
At iterate 22 f = 6055.3 |proj g|= 1.8976
At iterate 23 f = 6055.2 |proj g|= 1.6393
At iterate 24 f = 6055 |proj g|= 1.629
At iterate 25 f = 6054.6 |proj g|= 1.8972
At iterate 26 f = 6054.6 |proj g|= 1.5981
At iterate 27 f = 6054.5 |proj g|= 1.5928
At iterate 28 f = 6054.4 |proj g|= 1.5774
At iterate 29 f = 6054.2 |proj g|= 1.8957
At iterate 30 f = 6054.1 |proj g|= 1.8949
At iterate 31 f = 6052.6 |proj g|= 1.2987
At iterate 32 f = 6052.3 |proj g|= 1.2399
At iterate 33 f = 6051.9 |proj g|= 0.96675
At iterate 34 f = 6051.8 |proj g|= 0.96791
At iterate 35 f = 6051.8 |proj g|= 0.94212
At iterate 36 f = 6051.8 |proj g|= 1.8852
At iterate 37 f = 6051.8 |proj g|= 0.71337
At iterate 38 f = 6051.8 |proj g|= 0.71134
At iterate 39 f = 6051.7 |proj g|= 0.3276
At iterate 40 f = 6051.7 |proj g|= 0.19604
At iterate 41 f = 6051.7 |proj g|= 0.14445
At iterate 42 f = 6051.7 |proj g|= 0.13129
At iterate 43 f = 6051.7 |proj g|= 0.11073
At iterate 44 f = 6051.7 |proj g|= 0.65001
At iterate 45 f = 6051.7 |proj g|= 0.50626
At iterate 46 f = 6051.7 |proj g|= 0.21399
At iterate 47 f = 6051.7 |proj g|= 0.1184
At iterate 48 f = 6051.7 |proj g|= 0.11665
At iterate 49 f = 6051.7 |proj g|= 0.076832
At iterate 50 f = 6051.7 |proj g|= 0.11684
At iterate 51 f = 6051.7 |proj g|= 0.1167
At iterate 52 f = 6051.7 |proj g|= 0.09314
At iterate 53 f = 6051.7 |proj g|= 0.10182
At iterate 54 f = 6051.7 |proj g|= 0.050777
At iterate 55 f = 6051.7 |proj g|= 0.025595
At iterate 56 f = 6051.7 |proj g|= 0.11416
iterations 56
function evaluations 73
segments explored during Cauchy searches 68
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.114159
final function value 6051.72
F = 6051.72
final value 6051.720361
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6195.1 |proj g|= 1.8444
At iterate 1 f = 6190.4 |proj g|= 1.567
At iterate 2 f = 6180.2 |proj g|= 1.8736
At iterate 3 f = 6172.8 |proj g|= 1.8258
ys=-8.212e-01 -gs= 6.952e+00, BFGS update SKIPPED
At iterate 4 f = 6153.1 |proj g|= 1.7389
ys=-2.505e+00 -gs= 1.744e+01, BFGS update SKIPPED
At iterate 5 f = 6141.1 |proj g|= 1.694
At iterate 6 f = 6130.4 |proj g|= 1.9376
At iterate 7 f = 6113 |proj g|= 1.8015
At iterate 8 f = 6107.7 |proj g|= 1.8026
At iterate 9 f = 6085.2 |proj g|= 1.6781
At iterate 10 f = 6084 |proj g|= 1.8835
At iterate 11 f = 6082.6 |proj g|= 1.8806
At iterate 12 f = 6072.2 |proj g|= 1.8738
At iterate 13 f = 6070.1 |proj g|= 1.4886
At iterate 14 f = 6065.8 |proj g|= 1.381
At iterate 15 f = 6064.4 |proj g|= 1.3835
At iterate 16 f = 6064 |proj g|= 1.8899
At iterate 17 f = 6062.2 |proj g|= 1.8814
At iterate 18 f = 6061.5 |proj g|= 1.421
At iterate 19 f = 6059.1 |proj g|= 1.4174
At iterate 20 f = 6058.3 |proj g|= 1.4164
At iterate 21 f = 6056.7 |proj g|= 1.4007
At iterate 22 f = 6055.7 |proj g|= 1.901
At iterate 23 f = 6055.3 |proj g|= 1.4749
At iterate 24 f = 6054.9 |proj g|= 1.4674
At iterate 25 f = 6054.7 |proj g|= 1.8948
At iterate 26 f = 6054.5 |proj g|= 1.8955
At iterate 27 f = 6054.4 |proj g|= 1.8936
At iterate 28 f = 6053.8 |proj g|= 1.1737
At iterate 29 f = 6053.3 |proj g|= 1.8945
At iterate 30 f = 6053.2 |proj g|= 1.8917
At iterate 31 f = 6053.1 |proj g|= 1.353
At iterate 32 f = 6052.9 |proj g|= 1.4434
At iterate 33 f = 6052.8 |proj g|= 1.8925
At iterate 34 f = 6052.7 |proj g|= 1.441
At iterate 35 f = 6052.4 |proj g|= 1.4288
At iterate 36 f = 6052.3 |proj g|= 1.0089
At iterate 37 f = 6052.2 |proj g|= 1.8889
At iterate 38 f = 6052.1 |proj g|= 1.8888
At iterate 39 f = 6052.1 |proj g|= 1.6079
At iterate 40 f = 6052 |proj g|= 1.2944
At iterate 41 f = 6052 |proj g|= 1.1205
At iterate 42 f = 6052 |proj g|= 1.8863
At iterate 43 f = 6052 |proj g|= 1.429
At iterate 44 f = 6051.9 |proj g|= 0.92281
At iterate 45 f = 6051.9 |proj g|= 1.3317
At iterate 46 f = 6051.9 |proj g|= 1.3044
At iterate 47 f = 6051.9 |proj g|= 0.78826
At iterate 48 f = 6051.9 |proj g|= 0.80595
At iterate 49 f = 6051.9 |proj g|= 0.91726
At iterate 50 f = 6051.8 |proj g|= 1.8835
At iterate 51 f = 6051.7 |proj g|= 1.8839
At iterate 52 f = 6051.7 |proj g|= 0.62815
At iterate 53 f = 6051.7 |proj g|= 0.59062
At iterate 54 f = 6051.7 |proj g|= 0.09207
At iterate 55 f = 6051.7 |proj g|= 0.069011
At iterate 56 f = 6051.7 |proj g|= 0.016776
At iterate 57 f = 6051.7 |proj g|= 0.007937
iterations 57
function evaluations 64
segments explored during Cauchy searches 65
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.00793695
final function value 6051.72
F = 6051.72
final value 6051.720335
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6196.1 |proj g|= 1.9441
At iterate 1 f = 6158.6 |proj g|= 1.8479
At iterate 2 f = 6131.1 |proj g|= 1.8288
ys=-4.125e+00 -gs= 2.537e+01, BFGS update SKIPPED
At iterate 3 f = 6093.3 |proj g|= 1.7964
At iterate 4 f = 6091.6 |proj g|= 1.8966
At iterate 5 f = 6089.7 |proj g|= 1.7487
At iterate 6 f = 6078.4 |proj g|= 1.7556
At iterate 7 f = 6069.5 |proj g|= 1.7184
At iterate 8 f = 6066.2 |proj g|= 1.6136
At iterate 9 f = 6065.3 |proj g|= 1.9357
At iterate 10 f = 6063.3 |proj g|= 1.9207
At iterate 11 f = 6061.3 |proj g|= 1.909
At iterate 12 f = 6059.9 |proj g|= 1.8995
At iterate 13 f = 6057.7 |proj g|= 1.7114
At iterate 14 f = 6057.6 |proj g|= 1.7106
At iterate 15 f = 6057.1 |proj g|= 1.9015
At iterate 16 f = 6056.4 |proj g|= 1.9085
At iterate 17 f = 6056 |proj g|= 1.8998
At iterate 18 f = 6055.9 |proj g|= 1.6836
At iterate 19 f = 6055.7 |proj g|= 1.6795
At iterate 20 f = 6055.6 |proj g|= 1.8978
At iterate 21 f = 6055.5 |proj g|= 1.8959
At iterate 22 f = 6055.3 |proj g|= 1.6477
At iterate 23 f = 6055 |proj g|= 1.6138
At iterate 24 f = 6054.5 |proj g|= 1.5619
At iterate 25 f = 6054.1 |proj g|= 1.8977
At iterate 26 f = 6054 |proj g|= 1.8975
At iterate 27 f = 6053.4 |proj g|= 1.4395
At iterate 28 f = 6053 |proj g|= 1.889
At iterate 29 f = 6052.3 |proj g|= 1.8807
At iterate 30 f = 6052 |proj g|= 0.72993
At iterate 31 f = 6051.9 |proj g|= 0.84726
At iterate 32 f = 6051.8 |proj g|= 1.8843
At iterate 33 f = 6051.8 |proj g|= 0.47819
At iterate 34 f = 6051.8 |proj g|= 1.8846
At iterate 35 f = 6051.7 |proj g|= 0.70317
At iterate 36 f = 6051.7 |proj g|= 0.58719
At iterate 37 f = 6051.7 |proj g|= 0.53426
At iterate 38 f = 6051.7 |proj g|= 0.28181
At iterate 39 f = 6051.7 |proj g|= 0.27471
At iterate 40 f = 6051.7 |proj g|= 0.71023
At iterate 41 f = 6051.7 |proj g|= 0.7442
At iterate 42 f = 6051.7 |proj g|= 0.17579
At iterate 43 f = 6051.7 |proj g|= 0.12481
At iterate 44 f = 6051.7 |proj g|= 0.11666
At iterate 45 f = 6051.7 |proj g|= 0.80637
At iterate 46 f = 6051.7 |proj g|= 0.047588
At iterate 47 f = 6051.7 |proj g|= 0.041069
At iterate 48 f = 6051.7 |proj g|= 0.1168
At iterate 49 f = 6051.7 |proj g|= 0.11675
At iterate 50 f = 6051.7 |proj g|= 0.043049
At iterate 51 f = 6051.7 |proj g|= 0.042103
At iterate 52 f = 6051.7 |proj g|= 0.065384
At iterate 53 f = 6051.7 |proj g|= 0.11522
At iterate 54 f = 6051.7 |proj g|= 0.044015
iterations 54
function evaluations 68
segments explored during Cauchy searches 64
BFGS updates skipped 1
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0440146
final function value 6051.72
F = 6051.72
final value 6051.720388
converged
* The 4 best values (multistart) obtained are:
6051.72 6051.72 6051.72 6051.72
* The model corresponding to the best one (6051.72) is stored.
[1] 10.81038
[1] 369550.6
[1] 8.353776
[1] 5.220835
The problem with doing constraint first is that you end up with a non-hypercube shaped input space (corners have been knocked off by constraint), which is not ideal. We might therefore want different sensitivity measures for pre- and post-constrained ensemble.
[1] 8.353776
[1] 8.353615
[1] 5.220835
[1] 5.220903
It looks as though bwl_io is very influential across a number of variables, even though it didn’t appear that interesting in the parginal plots. Why is that?
A closer look at bwl_io now that the impact of f0_io has been removed shows a large number of “zero” NPP when bwl_io is at low values, which could well be the source of apparent sensitivity.
Take only higher values of bwl_io for the next round of constraints.
Define constraint level 1a as those members that run, and have F0_io <0.9 & b_wl_io > 0.15 (normalised).
# Plot the regular km emulator.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -5137.861
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 5137.9 |proj g|= 1.8356
At iterate 1 f = 5135.4 |proj g|= 1.7202
At iterate 2 f = 5131.9 |proj g|= 1.6
At iterate 3 f = 5131.2 |proj g|= 1.6826
At iterate 4 f = 5122.2 |proj g|= 1.6838
At iterate 5 f = 5119.7 |proj g|= 1.6697
At iterate 6 f = 5110.5 |proj g|= 1.6667
At iterate 7 f = 5105.6 |proj g|= 1.3939
At iterate 8 f = 5098.3 |proj g|= 1.4145
At iterate 9 f = 5093.3 |proj g|= 1.592
At iterate 10 f = 5091.8 |proj g|= 1.5626
At iterate 11 f = 5089.7 |proj g|= 1.4986
At iterate 12 f = 5088.5 |proj g|= 1.3214
At iterate 13 f = 5085.1 |proj g|= 1.501
At iterate 14 f = 5083.9 |proj g|= 1.4621
At iterate 15 f = 5082.1 |proj g|= 1.4517
At iterate 16 f = 5081.1 |proj g|= 1.494
At iterate 17 f = 5080.6 |proj g|= 0.95799
At iterate 18 f = 5080.4 |proj g|= 1.3623
At iterate 19 f = 5080 |proj g|= 0.56793
At iterate 20 f = 5079.9 |proj g|= 1.4786
At iterate 21 f = 5079.8 |proj g|= 0.46896
At iterate 22 f = 5079.8 |proj g|= 0.28689
At iterate 23 f = 5079.8 |proj g|= 0.57427
At iterate 24 f = 5079.7 |proj g|= 0.57268
At iterate 25 f = 5079.7 |proj g|= 0.36422
At iterate 26 f = 5079.7 |proj g|= 0.16654
At iterate 27 f = 5079.7 |proj g|= 0.22658
At iterate 28 f = 5079.7 |proj g|= 0.27311
At iterate 29 f = 5079.7 |proj g|= 0.21159
At iterate 30 f = 5079.7 |proj g|= 0.52762
At iterate 31 f = 5079.7 |proj g|= 0.2638
At iterate 32 f = 5079.7 |proj g|= 0.15332
At iterate 33 f = 5079.7 |proj g|= 0.16807
At iterate 34 f = 5079.7 |proj g|= 0.2156
At iterate 35 f = 5079.7 |proj g|= 0.093168
At iterate 36 f = 5079.7 |proj g|= 0.069853
At iterate 37 f = 5079.7 |proj g|= 0.11439
At iterate 38 f = 5079.7 |proj g|= 0.18854
At iterate 39 f = 5079.7 |proj g|= 0.14403
At iterate 40 f = 5079.6 |proj g|= 0.075684
At iterate 41 f = 5079.6 |proj g|= 0.045529
At iterate 42 f = 5079.6 |proj g|= 0.029163
At iterate 43 f = 5079.6 |proj g|= 0.028335
At iterate 44 f = 5079.6 |proj g|= 0.081292
At iterate 45 f = 5079.6 |proj g|= 0.03851
At iterate 46 f = 5079.6 |proj g|= 0.011842
At iterate 47 f = 5079.6 |proj g|= 0.010565
iterations 47
function evaluations 54
segments explored during Cauchy searches 60
BFGS updates skipped 0
active bounds at final generalized Cauchy point 20
norm of the final projected gradient 0.0105653
final function value 5079.65
F = 5079.65
final value 5079.648744
converged
It looks as though both the absolute value and the change over time are controlled by the same parameters.
The emulators appear to be at least capturing the broad response for all of the output variables.
First, plot the straight kriging emulators
[1] 5.433077
[1] 5.142409
[1] 5.2209
[1] 4.193005
[1] 4.219396
[1] 4.531841
[1] 5.629425
[1] 8.565018
[1] 5.224583
[1] 9.223968
[1] 7.089
[1] 7.78196
[1] 6.900971
[1] 6.929568
Next, plot the twostep glmnet/km emulators
We use the leave-one-out Mean Absolute Error, expressed as a percentage of the range of the output across the ensemble. We find that the twostep emulatorisn’t significantly more accurate, and is indeed less accurate for tree fraction.
[1] 5.410355
[1] 4.733929
[1] 5.215455
[1] 4.173501
[1] 4.03861
[1] 4.597044
[1] 5.501413
[1] 8.563912
[1] 5.213222
[1] 11.36184
[1] 7.260354
[1] 7.707121
[1] 6.725096
[1] 6.865819
To use History Matching, we need to specify targets for various model outputs. We treat these targets as “observations” with an uncertainty where a model run marked as “implausible” (beyond 3sd) matches the hard boundaries previously identified by A. Wiltshire as being desirable.
Choose the centre of the (implied) uniform distribution. cs_gb.target = (3000 - 750) / 2 = 1125 cv.target = (800 - 300) / 2 = 250 npp_n_gb.target = (80 - 35) / 2 = 22.5
nbp.target = 0 (gpp.target = 75) (runoff.target = 1)
(to do: visualise implausibility of design and loo emulated values of design as histogram)
Plot the original ensemble, plus the stuff that matches the constraints. Include anomalies. “Level 2” constraints are broad constraints produced by AW “Level 3” constraints are tighter constraints from citeable sources produced by ER.
We find that adding a constraint to level 2 constrains the behaviour of the soil carbon pool, but not the vegetation carbon pool. Probably not a surprise, given the Soil Carbon pool (and changes) are much larger than the vegetation carbon pool.
We can calculated the atmospheric growth as Fossil fule emissions (GCB) - Ocean uptake (GCB) - net land flux (JULES NBP)
Conversely, we can calculate the land uptake as NBP = FF - ocean - ag (Do we have AG in JULES?)
These two methods give very similar outcomes, bar an outlier.
We plot the JULES total land carbon change from 1850, together with “observations” from the Global Carbon Budget (Friedligstein et al, 2020). The “observations are the cumulative sum of the ‘land sink’ column on the”historical observations tab of the data from the paper. It states “The land sink is the average of several dynamic global vegetation models that reproduce the observed mean total land sink of the 1990s.”
I think to make a fair comparison, we might need to take the diff of the anomaly (or just the timeseries), as the above represents the total change since 1850, and I think all the “land sink” stuff looks at the year-to-year changes.
What do these parts of parameter space tell us about what we need to do to have the model reproduce reality?
The function addNroyDesignPoints builds an emulator for each model output in Y. It compares the output of each emulator at a number of candidate desin points, and chooses a space-filling set of them that that are Not Ruled Out Yet (statistically close to the observation at Y_target).
The function write_jules_design here simply takes the points calculated by addNroyDesignPoints and writes them to configuration files.
Check that the augmented design produces what we expect. New ensemble members should be somewhat constrained within the boundaries of the original design, if the comparison to data offers any constraint.
Visualise the predicted outputs at the NROY points of the old design, and the new suggested design.
Interestingly, these aren’t perfectly within the original hard boundaries set by Andy, even though I’ve set those boundaries to be the +- 4 standard deviation threholds in the History Match. I suggest this is because there is model discrepancy, and that there is considerable wriggle room induced from emulator uncertainty.
In particular, it appears that vegetation carbon is difficult to keep high, and that many NROY proposed members have a fairly low vegetation carbon. This might need a discrepancy term, or adjusting in some other way. It certainly needs exploring, and a OAAT plot might give clues as to the parameters to choose.
How good are the four emulators that we’ve built? Are there biases? (there’s no real evidence of this)
It’s obviously hard to maintain a high vegetation carbon in particular. What parameter values might you choose to do this, and what might be the trade-offs you have to make?
First, what parameters affect vegetation carbon and how? How sure are we about that?
What might be the trade-offs for a high (or accurate) vegetation carbon? are they acceptable? Plot the oaat sensitivity of the other 3 outputs we’re calibrating on.
What do the emulators make of the design points which actually make Andy’s “hard boundary” criteria? If we leave them out, do they still place the output within the hard boundaries?
Do a leave-one-out cross validation of points inside the hard boundaries using the wave1 fits.
We see in the leave-one-out analysis that the emulator is consistently under-predicting the vegetation carbon (though the uncertainty estimate often covers the actual value).This suggests (1) that there isn’t really a huge problem with a model discrepancy (or at least that isn’t the only problem), and (2) the history matching is working as it should, and taking into account a not-great emulator.
Is a two-step emulator any better at emulating those crucial points which fall within aw’s hard boundaries? First, create a list of emulator fits.
It doesn’t appear that the two-step emulator (here plotted in red) is doing any better than the regular emulator.
Two other things I can think of to check: 1) how about using “multistart” to choose different starting conditions for optimising the emulators and 2) Using a flat prior for the mean function rather than a linear prior.
Input from Eddy Robertson:
“We can assume that the majority of vegetation carbon is stored in tree trunks so, the carbon density of trees is approximately (cVeg_lnd_mean / treeFrac_lnd_mean). Although I suppose it’s possible that in the PPE the shrubs have become tree-like, so it might be interesting to plot cVeg_lnd_mean versus (shrubFrac_lnd_mean + treeFrac_lnd_mean) as well. The first question is whether cVeg_lnd_mean increases linearly with treeFrac_lnd_mean. Then if this turned out to be interesting, we could produce a treeCVeg_lnd_mean output.”
Is the range of carbon densities narrower in the “hard boundary” set?
For example, what effect has the AW “hard boundary” constraint produced on ranges of the other (marginal) outputs? Which constraint is doing the majority of the work?
A function to find the proportion of output space that is removed (or retained) when applying a constraint.
All constraints together take the normalized range of all the outputs to a little under a half of the original range, on average. Of course, four of the inputs are constrained here, which has quite a large impact.
We can see what impacts the constraints have individually
The NBP constraint has no marginal impact
Vegetation carbon (cVeg) has the strongest impact on everything else. This is no suprise, given that so few of the ensemble members seem to simulate carbon vegetation high enough.